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Abstract 

Generalizations of Bell polynomials, Bell numbers, and Stirling num- 
bers of the second kind have been introduced and their generating func- 
tions were evaluated. 
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1 Fractional Poisson Probability Distribution 

Bell polynomials, Bell numbers [T] and Stirling numbers of the second kind [5]- 
[4] are related to the well known Poisson probability distribution. Recently, the 
fractional Poisson probability distribution has been developed [5], and some of 
its applications have been implemented in [6] . The fractional Poisson probability 
distribution is a natural model which captures long memory impact on the 
counting process. 

Here we present and explore generating functions for new polynomials and 
numbers related to the fractional Poisson probability distribution. Those new 
polynomials and numbers are: fractional Bell polynomials, fractional Bell num- 
bers and fractional Stirling numbers of the second kind (see, [B]). 

The fractional Poisson probability distribution P^(n,t) of arriving n items 
(n = 0, 1, 2, ...) by time t is given by [5] 

, s {vt^Y (k + n)\ (-vt») k 
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where the parameter v has physical dimension [is] = sec _A1 and the gamma 

oo 

function has the familiar representation r(/i) = J dte^t^ 1 , Rc/i > 0. 

o 

At /i = 1 Eq.([T]) is transformed into the well known equation for the standard 
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Poisson probability distribution with substitution v — > u, where v is the rate of 
arrivals of the standard Poisson process with physical dimension v = sec -1 . 

2 Fractional Bell polynomials and Bell numbers 

Based on the fractional Poisson probability distribution Eq. (JTJ) , we introduce a 
new generalization of Bell polynomials 



where the parameter fi is < fj, < 1. We will call B^(x,m) as fractional 
Bell polynomials. Polynomials B^(x,m) are related to the well-known Bell 

oo 

polynomials [T] B(x,m) by B fl (x, m)| M=1 = B(x,jn) = e~ x ^ n m ^-j. From 

Eq.Q we come to the formula for numbers B^(m), which we call the fractional 
Bell numbers 



n m ^ (k + n)\ (-l) fe 

n=0 fc=0 u v ' " 

Now we focus on the general definitions given by Eqs.© and (J3j) to find 
the generating functions of polynomials B^(x,m) and numbers B^(m). Let us 
introduce the generating function F^(s,x) of polynomials B^{x,ni) as 

F ll (a,x) = J2— l B »(?>™)- ( 4 ) 

To find an explicit equation for F fJi (s,x), we substitute Eq.Q into Eq.Q 
and evaluate the sum over m. As a result we have 

F ll (s,x) = E fi (x(e s -l)), (5) 
where E^z) is the Mittag-LefHer function given by its power series [7] 

oo m 

If we put x = 1 in Eq.([5]), then we immediately come to the generating 
function B^{s) of the fractional Bell numbers B^{ni) 



W = E^H-^(e S -l). (7) 



m=0 
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In the case of /i = 1, Eq.© turns into the equation for the generating 
function of the Bell polynomials Fi(s,x) — exp{ir(e s — 1)}, while Eq. dTJ) reads 

^i( s ) = S ^—\Bi(m) — exp(e s — 1), and we come to the equation for the Bell 

m=0 

numbers generating function. 



3 Fractional Stirling numbers of the second kind 

Now we introduce the fractional generalization of the Stirling numbers of the 
second kind S^{m,l) by means of equation 



B^(x, m) = ^2 Sn(m, l)a 



(8) 



where B^(x,m) is a fractional generalization of Bell polynomials given by 
Eq.© and the parameter /it is < fj, < 1. At /i — 1, Eq.© defines integers 
S(m,l) = S^(m, /)|ft=i) which are called Stirling numbers of the second kind. 
At x = 1, when fractional Bell polynomials B^{x,m) become fractional Bell 
numbers, B^(m) = B ^(x , m)\ x= i , Eq.© gives us the equation to express frac- 
tional Bell numbers in terms of fractional Stirling numbers of the second kind 

m 

B n{m) = 2 Sn(m,l). 
1=0 

To find Sfj,(m,l) we transform the right-hand side of Eq.© as follows 



s.fe".) = Ero5TT)E<-i)'-" 

1=0 yr ' n=0 
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where the notation (^) = ^mz^\ has been introduced. By comparing Eq.© 
and Eq.© we conclude that the fractional Stirling numbers S^(m,l) are given 

by 



(10) 



S M (rn, 0) = 5 m fi, 5^(m,/) = 0, l = m+l, m + 2, .... 

As an example, Table 1 presents a few fractional Stirling numbers of the 
second kind. 
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r(p+i) 



r(/H-i) 



r(M+i) 



~r 

FES 



r(p+i) 



r(2M+i) 



r(2^+i) 



~TT~ 



r(2H-i) 



~3TT 

r(2^+i) 



r(2 M +i) 



r(3H-i) 



T(3m+1) 



T5TT 

T(3m+1) 



51TT" 

r(3M+l) 



"2T" 



r(4M+l) 
240 



r(4M+i) 

1560 

r(4 M +i) 



120 

r(5M+i) 

1801) 

t(5m+i) 



T(6m+1) 



Table 1. Fractional Stirling numbers of the second kind S^m, I) (0 < fj, < 1) 
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To find a generating function of the fractional Stirling numbers S^m, I), let's 
substitute B^(x,m) from Eq.® into the definition given by Eq.Q. Hence, we 
have 

00 m / m \ 00/00 m \ 

m=0 ' \/=o / ;=o \m=i '/ 

On the other hand, from Eq.©, we have F^(s,x) — w pnr ^ ■ Upon 

z=o 

comparing this equation and Eq.(|lip we introduce the generating function {? M (s, Z) 
0,(8, 1) = £ S„(m, = ffe^ * = °> L 2 ' •••• ( i2 ) 

m— / 

In addition, we introduce the generating function J-^(s,t) defined by 

00 m m 1 00 j, s w 

t) = EE s -M^-r = E Y^-yV = - !))• (13) 

m=0 i=0 ' (=0 VM ' 

As a special case /1 = 1, equations (|12|) and (fT3|) include the well-know 
generating function equations for the standard Stirling numbers of the second 
kind S(m,l) (for example, see Eqs.(2.17) and (2.18) in Ref.g]). 
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